Draft version January 18, 2012 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



DARK MATTER DECAY AND ANNIHILATION IN THE LOCAL UNIVERSE: CLUES FROM FERMI 



A. J. CUESTA 



(N 

o 



X 
6 



> 

m 

o 
o 



% 



T. E. Jeltema^, F. Zandanel^ S. Profumo'*'^, F. Prada^"', G. Yepes^', A. Klypin^, Y. Hoffman*^ 
S. Gottlober'^ J. Primack^, M. A. Sanchez-Conde"'-" and C. Pfrommer^^ 

Draft version January 18, 2012 



ABSTRACT 

We present all-sky simulated Fermi maps of 7-rays from dark matter decay and annihilation in the 
Local Universe. The dark matter distribution is obtained from a constrained cosmological simulation 
of the neighboring large-scale structure provided by the CLUES projeco The dark matter fields of 
density and density squared are then taken as an input for the Fermi observation simulation tool to 
predict the 7-ray photon counts that Ferm,i would detect in 5 years of all-sky survey for given dark 
matter models. Signal-to-noise sky maps have also been obtained by adopting the current Galactic and 
isotropic diffuse background models released by the Ferm,i collaboration. We point out the possibility 
for Fermi to detect a dark matter 7-ray signal in local extragalactic structures. In particular, we 
conclude here that Fermi observations of nearby clusters (e.g. Virgo and Coma) and filaments are 
expected to give stronger constraints on decaying dark matter compared to previous studies. As an 
example, we find a significant signal-to-noise ratio in dark matter models with a decay rate fitting the 
positron excess as measured by PAMELA. This is the first time that dark matter filaments are shown 
to be promising targets for indirect detection of dark matter. On the other hand, the prospects for 
detectability of annihilating dark matter in local extragalactic structures are less optimistic even with 
extreme cross-sections. We make the dark matter density and density squared maps available online 
at http://www. clues-project . org/articles/darkmattermaps .html. 

Subject headings: astroparticle physics — dark matter — large-scale structure of Universe — gamma 
rays: diffuse background — methods: numerical 



1. INTRODUCTION 

A large amount of astrophysical evidence suggests that 
most of the Universe's matter content is in the form of 
cold dark matter (DM). However, the nature of DM is 
still one of the most important open questions in modern 
physics. Many different candidates have be en proposed 
as DM constituents (see iBertone et al.|[2005l for a review 
on candidates, and experimental searches), yet for the 
time being there is no evidence in favor of any model. 
One of the most studied scenarios is that of weakly inter- 
acting massive particles (WIMPs) , where 7-rays are gen- 
erated as secondary products of WIMP decay or annihila- 
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tion (e.g. lBertondl2010l ). Therefore, 7-ray observations, 
being a complementary approach to direct searches, are 
a powerful tool to study the nature of DM. 

At present, the Imaging Atmospheric Cerenkov Tele- 
scopes (lACTs such as MAGIC, HESS and VERITAS) 
together with the recently launched Fermi satellite of- 
fer great tools to search for the 7-ray emission due 
to DM decay or annihilation in the MeV-TeV energy 
range. The main instrument on board Ferm,i is the 
Large Area Telescope (LAT), which is designed to ex- 
plore the en tire 7-ray sky in th e 20 MeV-300 GeV en- 
ergy range (jAtwood et aLll2009f l. The Fermi-LAT col- 
laboration already reported some of their results on 
DM searches, with no detection of DM 7-ray e mis- 
sion in dwarf spheroidal galaxie s ([AbdoIeranilQlOc), 
clusters (lA ckcrma nn et aiTl2 010'l. or spectral features 
(|Abdo et al. 2010a. lAbdoet al. 2010ta' ) lACTs did not 
succeed in DM detection either fe g lAharonian et al. 



I2008al lAliu et al.l[2009l . lAleksic et al]l2010LIAcciari et al 
2010). These studies are focused mainly on DM annihi- 
lation. 

Despite these negative results, the recent detec- 
tion of an excess of high-energy (10-100 GeV) 
positrons over the standard expectation from galac- 
tic cosmic-ray mode ls by the PAMELA experiment 
(|Adriani et al.l [2009b '). triggered an interest in the pos- 
sibility that these positrons originate from DM in the 
Milky Way. Interestingly, if DM annihilation or de- 
cay is to consistently explain the cosmic ray elec- 
tron/positron data recently produced by th e PAMELA, 

Ferrn i-hKI a nd HES S Collaboration s (lAbdo et al.l 

20091 Aharon ian et al.l [2008 b. Aharoni an et all I2009I . 



Adriani et al.ll2009al . lAdriam et al.ll20ig ). one is generi- 
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cally forced to consider dominantly leptonic final states 
to avoid over-producing anti-protons and to generate 
enough high-energy positrons; and a cross-section that 
is enhanced over its standard value by ^ 10'^. In 
turn, this is an interesting scenario for 7-ray searches, 
since hard final state leptons yiel d an unmistakable hard 
bremsstrahlung 7-ray spectrum (jPinzke et al.ll20d9l ). 

The aim of this letter is to study DM induced 7-ray 
emission using a constrained Local Universe cosmologi- 
cal simulation, as would be observed by the Fermi satel- 
lite in an all-sky 7-ray survey for given DM models. In 
Section [21 we describe the cosmological simulation used 
to infer the DM distribution in the Local Universe. In 
Section |3l the Fermi observation simulations and the as- 
sumed particle physics models are described. In Section|4] 
we present the main results from the 7-ray all-sky maps. 
We finally discuss the main conclusions of our work in 
Section [5] 

2. CONSTRAINED SIMULATIONS OF THE LOCAL 
UNIVERSE 

In order to get a detailed description of the DM den- 
sity distribution in the Local Universe, we use a high- 
resolution cosmological simulation box from the CLUES 
Project. This simulation set provides a realistic local 
density fiel d which is consist ent with the ACDM cos- 
mology fsee lYepes eraIll2009l and lGottl6eber et al.ll2010l 
for more details). As we want to study 7-rays from 
large structures in the Local Universe such as nearby 
galaxy clusters, we choose the Boxl60CR simulation. 
This is a constrained realization with 1024^ particles in 
a cube of 160/i^^Mpc on a side which was run using 
the MPLART cosmologica l code (|Kravtsov et al.lll997L 
iGottloeber fc KlvpinI I2008D . The initial conditions are 
set assuming WMAP3 cosmology (with f2,„ = 0.24, 
flA = 0.76, nt = 0.042, h = 0.73, erg = 0.75, and 
n = 0.95) and implements the constraints from the ob- 
served density field so that it reproduces the observed 
matter distributio n in the Local Universe on large scales 
at re dshift z = (jHoffman fc Ribak|[T99l iKlvpin et all 
120031) . The massive clusters such as Virgo, Coma and 
Perseus, together with the Great Attractor, are well re- 
produced. However, the final positions of these objects 
are not exactly at their observed positions, with a typical 
error around 5/i~^Mpc. 

This cosmological simulation allows us to produce all- 
sky maps of the local DM density and density squared, 
which are proportional to 7-ray emission due to particle 
DM decay and annihilation, respectively (to include more 
distant stru ctures, one could u se a box replication tech- 
nique as in iZavala et al.ll2010D. We follow the method 
described in iKuhlen et al.l ( )2008D to compute these lu- 
minosities. The flux is proportional to ^^ mp/Andf for 
decay and ^^ mpPi/Airdf for annihilation, where i runs 
from 1 to the number of particles in each pixel, rup is 
the mass of the simulation particle, pi is the density as- 
sociated to the i-th particle, computed using the sphere 
which contains its 32 nearest neighbors (no smoothing 
kernel was used), and di is the distance to the observer. 
Only particles between 5/i~^Mpc and 80/i~^Mpc from 
the observer, which is placed at the right distance from 
the Virgo cluster, are taken into account. A proper de- 
scription of the density field in the innermost 5h~^Mpc, 
although affected by random density fiuctuations at this 



scale, would require a higher resolution simu lation to be 
resolved properly (as in lLibeskind et al.]|2010D . Thus this 
region empty of massive large structures is not considered 
here. We bin these fiuxes in a Cartesian grid with 3600 
and 1800 pixels of galactic longitude and latitude, re- 
spectively. This corresponds to an angular resolution of 
roughly 0.1 deg per pixel, reproducing the best angular 
resolution that Fermi has at its highest accessible energy 
range. 

Due to the finite resolution of the simulation, we can- 
not resolve the very inner center of DM halos. For this 
reason, we correct the fiux, if underestim ated, in every 
pixel where the centers of DM halos lie (see | Kuhlen et al.l 
I2008D . We assume a NFW profile ([Navarro et al.lll996|) 
for the inner density profile of these halos and we extrap- 
olate it up to the halo center. The scale radius r^ of these 
halos is calculated from the virial mass-concentration re- 
lation in Maccio et aLs (2008). Typical corrections do not 
exceed ^^25% and ~250% of their original value for de- 
cay and annihilation, respectively. We note that no boost 
factor due to DM substructures, or any other effect (such 
as adiabatic contraction from baryons or Sommerfeld en- 
hancement), is included in our analysis. 

The resulting all-sky maps are shown in Figure [U 
where known objects are highlighted. The images are 
done with HEALPi^fl. Angular projection used here is 
equirectangular (plate carree) . Maps are color-coded ac- 
cording to the logarithmic flux in each 0.1 deg pixel, mea- 
sured in GeV/c^ cm~^ kpc sr"-'^ for decay and GeV^/c'* 
cm~^ kpc sr~^ for annihilation. These maps are used as 
input for the Fermi-hAT observation simulations, which 
we describe in the following Section. 

3. FERMI SATELLITE OBSERVATION SIMULATIONS 

Using the full sky DM density and density squared 
maps, the simulated Fermi-LAT observations are pro- 
duced using the gtobssim routine, part of the Fermi 
Science Tools package (v9rl5p2), which incorporates the 
Fermz-LAT effective area and point spread function and 
their energy dependence. All simulations are run to gen- 
erate a 5-year observation in the default scanning mode 
and using the current release of the LAT instrument re- 
sponse functions (P6.V3_DIFFUSE). 

In the present study, we specifically adopt two exam- 
ples for the 7-ray spectrum from decay or annihilation of 
the DM particle, chosen to be representative of more gen- 
eral classes of DM models. The first model features a DM 
particle with a inass of 1.6 TeV yielding a pair oi fi~^fi~ 
([Bergstrom et al.ll2009( ). wh ich was shown to fi t accu - 
rately the PAMELA data in [Papucci fc Strumial ([2010D . 
In this case, 7-ray emission is produced directly in the 
final state radiation (FSR) as well as through inverse 
Compton scattering (IC) of the high energy e"*" and e~ 
produced off of cosmic microwave background photons. 
We include the expected contributi ons from both FSR 
and IC (see [Ackermann et al.[ [20T0I for details). Given 
that the energy loss time scales for high-energy elec- 
trons and positrons produced by muon decays are much 
shorter than the diffusion time scales in the structures 
we consider here, we neglect diffusion, and also calcu- 
late the emission of said electrons and positrons via IC 
up-scattering of CMB photons. This yields a significant 

^ |http : //healpix ■ jpl ■ nasa . gov] 
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low-energy component, extending al l the way up to ener- 
gies r elevant to the Fermi telescope (jProfumo fc Jelteiiial 
l2009f ). We also consider a second, more conventional 
model, inspired by what it is expected in e.g. super- 
symmetric models with a bino-like lightest supersym- 
metric particle: a 100 GeV neutralino yielding a quark- 
antiquark pair (of b flavor, for definiteness). The primary 
source of 7-rays here is the decay of neutral pions pro- 
duced in the bb hadronization chains. 

The Fermi simulations of the 7-ray signal from DM 
annihilation (both to bb and fj.'^fi^) were normalized to 
a DM flux in the Fermi energy range of 9 x 10~* photons 
cm~^ s~^ integrated over the full sky. This was chosen 
in order to obtain good statistics to compare between 
different extragalactic structures. For bb, this flux cor- 
responds to a cross-section of 10~^^cm^s~^. This cross- 
section value has already been excluded by current indi- 
rect searches with Fe rmi, PAMELA and HESS (see e.g. 
ICireUi et ahllMol and Fapucci fc Strumiall2010l ). Yet we 
decided to keep this extreme case for comparison pur- 
poses. For ^"""/i", the same total flux corresponds to an 
annihilation cross-section of 5.8 x lO^^'^cm'^s"^, which 
gives a good fit to the PAMELA posi tron excess. This 
value is only marginally excluded in lAckermann et aLl 
(120101) if cluster substructur es are considered, and also 
in lPapucci fc Strumial (|2010D if the Milky Way DM halo 
follows NEW. 

In the case of DM decay as e.g. in supersymmetry with 
very weak R-parity violation jnduced by a dimension- 6 
GUT-scale operator to both bb and /x^/x", we simulated 
a total Fermi flux over the full sky of 1.5 x 10~^ photons 
cm~^ s~^, which corresponds to a decay lifetime of t ~ 
10^^ s for bb and t ~ 3 x 10^^ s for ^J.'^fJ,^ ■ These hfetimes 
are not currently excluded by other 7-ray constraints and 
the latter case gives a good fit to^the PAM ELA measured 
positron fraction e xcess (e.g. 



Cirelli et al. 



Zhang et al 



2010) 



|P; 
II 



Papucci & Strumia 



2OIOI: IChen et al.ll20ia .Meade et al 



2010; 



20101: 



We also include in the simulations realistic treatments 
of both Galactic and isotropic diffuse backgrounds. In 
particular, the 7-ray emission from the Galaxy is quite 
variable across the sky, an important consideration when 
comparing the expected signals from known objects. Eor 
example, structures lying at low Galactic latitudes like 
the Great Attractor will have much higher 7-ray back- 
grounds than high latitude objects like the Virgo cluster. 
We simulated 5-year Fermi observations of the Galac- 
tic and isotropic diffuse backgrounds using the current 
background models released by the Fermi collaboratiorQ 
(gll_iem_v02 .f it and isotropic_iem_v02.txt, respec- 
tively). The output background map is then used to 
compute signal-to-noise (S/N) all-sky maps, detailed in 
the next Section. 

4. RESULTS 

In Figure [3] we show our main result: the S/N pre- 
diction for the extragalactic 7-ray emission in the 100 
MeV-10 GeV energy range from annihilation and decay 
of DM in the Local Universe, as it would be seen by 
the Fermi satellite after 5-years of observations. This 
is the first time that a constrained cosmological simula- 

^ http : //f ermi . g sf c .nasa . gov/ssc/data/access/lat/ | 
BackgroundModels . html 



tion is used to generate maps that are consistent with 
both the currently accepted cosmology and the observed 
Local Universe. These maps assume a particle rnass of 
100 GeV which annihilates or decays through the bb chan- 
nel. Maps for the DM model yielding /^"'"/i" are simi- 
lar but present lower S/N. Pixels are binned in squares 
of 1 deg which matches the Ferm,i-LAT PSF at around 
1 GeV as well as the typical angular size in the sky of 
nearby clusters. In order to make a quantitative analysis 
of DM detectability in large nearby structures, including 
cluster and filament regions, we computed in Table [5] the 
photon number counts and S/N from annihilation and 
decay in the 1 GeV-10 GeV energy range, for both the 
bb and iJ.'^fJ.^ channels. The S/N is defined as the signal 
over the square root of the signal plus the background 7- 
ray emission. This choice of energy range maximizes the 
S/N ratio as compared to the 100 MeV-10 GeV range, as 
the fermi-LAT sensitivity is significantly worse at lower 
energies. 

In the case of DM decay, we find that nearby clus- 
ters and filamentary regions could be detected for decay 
lifetimes longer than those currently ruled out by other 
7-ray constraints as discussed above. This shows that 
extragalactic structures are excellent targets to search 
for a signal or to place constraints on DM decay models, 
including those fitting the PAMELA positron data. It 
is important to note that according to Table [2] the most 
promising clusters for DM studies are high galactic lati- 
tude objects, like Virgo and Coma. The former, however, 
has not been used in the recent Fermi search for DM in 
clusters, due to the presence of M87, which is a powerful 
7-ray source, which makes the extraction of a signal or 
a limit much more complicated than for other clusters. 
Moreover, we find that S/N is not very sensitive to the 
area of the region under analysis, provided that the aper- 
ture radius is no more than few degrees, where the signal 
saturates and therefore the background noise makes the 
S/N decrease. 

We also highlight that, in the case of DM decay, the 
filamentary structure of the cosmic web constitutes an in- 
teresting target for DM searches. To our knowledge, this 
is the first time that filaments have been considered as 
targets for DM searches. In this case the 7-ray luminos- 
ity is just proportional to the enclosed mass, where as this 
is onl y approximately true for annihilation (Pinzk e et al.l 
I2009f ). This means that massive extragalactic objects of- 
fer the best chance for detection (see Table ^ . Large 
filaments of DM match and even exceed the values of 
S/N as compared to those in large clusters, although cau- 
tion should be taken regarding the exact orientation in 
the sky of these filaments as we find some variation in 
smaller volume constrained simulations. Superclusters 
such as the region marked in Figure [3] show even more 
significant values. Hence, these features of the Large- 
Scale Structure of the Universe may prove to be a very 
promising novel way to detect decaying DM with Fermi. 

5. DISCUSSION AND CONCLUSION 

In this letter we have presented simulated Fermi 
maps of extragalactic 7-rays coming from DM anni- 
hilation and decay in the Local Universe. The DM 
distribution is taken from a constrained cosmological 
simulation of the Local Universe Boxl60CR by the 
CLUES project, and it is available online at the follow- 
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Fig. 1. — DM distribution in the Local Universe constrained cosmological simulation Boxl60CR. These all-sky maps are Carte- 
sian projections in Galactic coordinates. Top panel shows the density distribution, whereas the bottom panel displays the distribu- 
tion of density squared. The maps are color-coded according to the logjg of the DM flux, and units are GeV/c^ cm"'' kpc sr~^ 
for decay map and GeV^/c* cm~^ kpc sr~^ for the annihilation map. Largcstructures reproduced by the simulation such as Virgo, 
Coma, and Perseus clusters, together with the Great Attractor, are labeled. High-resolution versions of these figures are available at 
Ihttp: //www. clues-project .org/articles/darkmattermaps .html 
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Fig. 2. — S/N all-sky maps from Fermi simulations for DM 7-rays in the energy range 100 MeV-10 GeV built from the Boxl60CR 
constrained simulation of the Local Universe. Results for DM decay (top) and annihilation (bottom) are shown for the bb channel model. 
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TABLE 1 

PHOTON COUNTS AND BACKGROUNDS IN Fermi SIMULATIONS. 



Object 




b b channel 
ann dec 


fi+fi 
ann 


channel 
dec 


background 


Coma 1 deg 
Coma 2 deg 
Coma 5 deg 






2.981 
2.592 

1.542 


(20) 
(31) 
(39) 


2.981 
4.061 
4.976 


(20) 

(52) 

(135) 


0.567 
0.371 
0.203 


(3) 
(4) 
(5) 


1.078 
1.331 

1.542 


(6) 
(15) 
(39) 


25 
112 
601 


Virgo 1 deg 
Virgo 2 deg 
Virgo 5 deg 






2.739 

2.858 
1.818 


(15) 
(28) 
(42) 


3.133 
5.371 
7.700 


(18) 

(61) 

(203) 


0.000 
0.471 
0.269 


(0) 
(4) 
(6) 


0.485 
1.444 
2.189 


(13) 
(51) 


15 

68 

492 


Perseus 1 deg 
Perseus 2 deg 
Perseus 5 deg 






1.144 
0.700 
0.386 


(12) 
(14) 
(20) 


5.642 
5.646 
4.542 


(74) 
(128) 
(245) 


0.493 
0.253 
0.116 


(5) 
(5) 
(6) 


1.841 
1.518 
1.187 


(20) 
(31) 

(62) 


98 

386 

2665 


GAttractor 1 deg 
GAttractor 2 deg 
GAttractor 5 deg 




0.280 
0.211 
0.136 


(S) 
(13) 
(23) 


2.686 
2.581 
2.157 


(80) 
(162) 
(367) 


0.105 
0.049 
0.041 


(3) 
(3) 
(7) 


0.935 
0.696 
0.538 


(27) 
(43) 
(91) 


807 

3777 

28572 


Filament 1, d = 


= 65Mpc/h 


0.224 


(14) 


4.515 


(292) 


0.128 


(8) 


1.348 


(85) 


3891 


Filament2, d = 


= 40Mpc//i 


0.752 


(67) 


9.317 


(871) 


0.191 


(17) 


2.589 


(233) 


7869 


t-'ilamentS, d = 


= 65Mpc//i 


0.351 


(84) 


4.862 


(1174) 


0.121 


(29) 


1.181 


(283) 


57127 


Filament4, d = 


= 55Mpc/;i 


0.576 


(91) 


8.380 


(1358) 


0.184 


(29) 


2.065 


(328) 


24904 


Superclusterl, 


d = 


45Mpc//i 


0.911 


(144) 


12.598 


(2066) 


0.254 


(40) 


3.334 


(531) 


24829 



Note. — The S/N ratio and number of photon counts (in brackets) in the 1 GeV— 10 GeV energy range 
for our different DM models. For cluster regions, three different radii are considered (1, 2, and 5 degrees). 
Filaments 1 to 4 represent elongated regions connected to these clusters which are potentially interesting 
due to their high S/N. Median distance of halos belonging to these filaments is indicated. Superclusterl is a 
collection of massive halos which accidentally lie along the line-of-sight. Background counts from the Galactic 
plus extragalactic diffuse in the same regions are also listed. Note that the annihilation to 6 6 case is shown 
for comparison purposes only. 



ing URL: http://www.clues-project.org/articles/ 
darkmattermaps.html. This distribution is then taken 
as an input for the Fermi observation simulation to ob- 
tain the all-sky distribution of 7-ray photon counts that 
would be measured by Fermi in a 5-year survey. Galac- 
tic and isotropic 7-ray diffuse backgrounds are also taken 
into account. This allows us to get S/N all-sky maps to 
estimate the possibility of detection of 7-rays from DM in 
such a survey. We adopted two different particle physics 
models: a DM particle annihilating (decaying) primarily 
to a bb final state, and a DM model that gives a good 
fit to the local positron fraction measured by PAMELA 
and the total electron spectrum measured by Fermi with 
annihilation (decay) to a fJ-'^ IJ.^ final state. 

In the case of DM decay, we find for these models that 
large clusters with high galactic latitudes, like Virgo and 
Coma, offer the best chance to be detected together with 
filamentary regions and large superclusters. Besides, we 
find that the S/N is not a strong function of the area of 
the analyzed region, which allows for considering large 
apertures without significant penalty in the results by 
masking away 7-ray point-sources. This is an important 
result as the Fermi collaboration and other authors have 
started to severely constrain models of annihilating DM 
Ackcrmann et al."2010l JAbdo et alll2010al . lAbdo et all 
2010b, Abdo et al. 201(3) while decaying DM has com- 
parative ly received consid erably smaller attention (how- 
ever, see lHiitsi et al.ll201Cl( ). Currently available 7-ray ob- 
servations provide less stringent constraints in this case, 
and from a theoretical standpoint, decaying DM is a 
generic prediction of many theories beyond the Standard 
Model of particle physics, for both neutralino DM and 
models that explain the PAMELA positron excess. 

We do not find any strong arguments in favor of a 
possible Fermi detection of extragalactic 7-rays induced 
by DM annihilation. However, we cannot exclude com- 






pletely this possibility since we are not considering possi- 
ble boost factors from DM substructures, adiabatic com- 
pression, and Sommerfeld effect that may significantly 
enhance the final 7-ray emission. Besides, the predicted 
signals will be enhanced due to recent determinations 
of the cosmological parameter ag suggesting a higher 
value than the one assumed here, which is consistent with 
WMAP3 cosmology 

These conclusions should be complemented by an anal- 
ysis of the galactic components, mainly from DM sub- 
halos like those hosting the DM-rich dwarf spheroidal 
galaxies arou nd the Milky Way. T his has been recently 
addressed in lAnderson et al.l ()2010l) based on the anal- 
ysis of a high-resolution sim ulation of a g a lactic DM 
halo, as in previous papers bv lKuhlen et al.1 (|2008h and 
iSpringel et al.1 (|2008D . We note that the presence of any 
galactic foregrounds not modeled here has a potential 
effect on the significance of our predictions, although 
only the Galactic Center and massive subhalos have been 
shown in these papers to be relevant. Moreover, the an- 
nihilation or decay of DM in Galactic subhalos will pro- 
duce 7-ray photons similar to those from DM in local 
extragalactic structures and thus if these happened to be 
coincident in the sky it would only enhance the signal. 
Nevertheless, a spectral confirmation of the potential DM 
signal is necessary to validate any claim of detection. On 
the other hand, these results on DM search will benefit 
from additional hints from the study of the angular power 
spectrum of the 7 -ray flu x (e.g. iFornasa et al.l 1 20091 
iHenslev et al.ll2009l see also lZavala et al.ll2010l for a sim- 
ilar approach to that presented here, extended to more 
distant contributions). An analysis of the anisotropics in 
the extragalactic diffuse radiation from DM using con- 
strained simulations will be presented elsewhere. 

Concluding, we find that Fermi will be able to place 
strong constraints on the nature of dark matter by 
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studying extragalactic objects, in particular for decay. 
The theoretical predictions from constrained simulations 
should provide the astroparticle community with the 
most interesting prospects for the detection of the 
elusive DM particle. 

We thank support of Spanish MICINNs Consolider- 
Ingenio 2010 Programme MULTIDARK CSD2009-00064 
and ASTROMADRID (S2009/ESP-146). A.J.C. thanks 
support from MEC Spanish grant FPU AP2005-1826. 



GY acknowledges support of MICINN research grants 
FPA2009-08958 and AYA2009-13875-C03-02. SP is sup- 
ported by NASA, DoE and NSF. YH is supported by Is- 
rael Science Foundation (13/08). SG acknowledges sup- 
port of DAAD through PPP program. BOX160GR sim- 
ulation has been performed at Leibniz Rechenzentrum 
Munich (LRZ). Fermi-L AT simulation tools were pro- 
vided by the Fermi-ljAT collaboration and the Fermi 
Science Support Center. 



REFERENCES 



Abdo, A. A., et al. 2009, Physical Review Letters, 102, 181101 

— . 2010a, Journal of Cosmology and Astro-Particle Physics, 4, 14 

— . 2010b, Physical Review Letters, 104, 091302 

— . 2010c, ApJ, 712, 147 

Acciari, V. A., et al. 2010, ApJ, 720, 1174 

Ackermann, M., et al. 2010, Journal of Cosmology and 

Astro-Particle Physics, 5, 25 
Adriani, O., et al. 2009a, Physical Review Letters, 102, 051101 
— . 2009b, Nature, 458, 607 
— . 2010, Physical Review Letters, 105, 121101 
Aharonian, P., et al. 2008a, Astroparticle Physics, 29, 55 
— . 2008b, Physical Review Letters, 101, 261104 
— . 2009, A&A, 508, 561 
Aleksic, J., et al. 2010, ApJ, 710, 634 
Aliu, E., et al. 2009, ApJ, 697, 1299 
Anderson, B., Kuhlen, M., Johnson, R., Madau, P., &; Dicmand, 

J. 2010, ApJ, 718, 899 
Atwood, W. B., et al. 2009, ApJ, 697, 1071 
Bergstrom, L., Edsjo, J., & Zaharijas, G. 2009, Physical Review 

Letters, 103, 031103 
Bertone, G. 2010, Cambridge University Press 
Bertone, G., Hooper, D., & Silk, J. 2005, Phys. Rep., 405, 279 
Chen, C, Mandal, S. K., & Takahashi, F. 2010, Journal of 

Cosmology and Astro-Particle Physics, 1, 23 
Cirelli, M., Panel, P., &; Serpico, P. D. 2010, Nuclear Physics B, 

840, 284 
Fornasa, M., Fieri, L., Bertone, C, & Branchini, E. 2009, 

Phys. Rev. D, 80, 023518 
Gottloeber, S., Hoffman, Y., & Yepes, G. 2010, High Performance 

Computing in Science and Engineering, Wagner, S., Steinmetz, 

M., Bode, A., &; Miiller, M.M. (Eds.), Springer- Verlag, 309 
Gottloeber, S., & Klypin, A. 2008, High Performance Computing 

in Science and Engineering, Wagner, S., Steinmetz, M., Bode, 

A., & Brehm, M. (Eds.), Springer- Verlag, 29 



Hensley, B. S., Siegal-Gaskins, J. M., & Pavlidou, V. 2009, 

arXiv:0912.1854 
Hoffman, Y., & Ribak, E. 1991, ApJ, 380, L5 
Hiitsi, G., Hektor, A., &; Raidal, M. 2010, Journal of Cosmology 

and Astro-Particle Physics, 7, 8 
Klypin, A., Hoffman, Y., Kravtsov, A. V., & Gottlober, S. 2003, 

ApJ, 596, 19 
Kravtsov, A. V., Klypin, A. A., & Khokhlov, A. M. 1997, ApJS, 

111, 73 
Kuhlen, M., Diemand, J., & Madau, P. 2008, ApJ, 686, 262 
Libeskind, N. I., Yepes, G., Knebe, A., Gottlober, S., Hoffman, 

Y., & Knollmann, S. R. 2010, MNRAS, 401, 1889 
Maccio, A. V., Dutton, A. A., & van den Bosch, F. C. 2008, 

MNRAS, 391, 1940 
Meade, P., Papucci, M., Strumia, A., & Volansky, T. 2010, 

Nuclear Physics B, 831, 178 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 

563 
Papucci, M., & Strumia, A. 2010, Journal of Cosmology and 

Astro-Particle Physics, 3, 14 
Pinzke, A., Pfrommer, C, & Bergstrom, L. 2009, Physical Review 

Letters, 103, 181302 
Profumo, S., &; Jeltema, T. E. 2009, Journal of Cosmology and 

Astro-Particle Physics, 7, 20 
Springel, V. et al. 2008, Nature, 456, 73 
Yepes, G., Martinez- Vaquero, L. A., Gottlober, S., & Hoffman, Y. 

2009, in AIPC Series, Vol. 1178, The CLUES project: 
Constrained Local UniversE Simulations, ed. C. Balazs & 
F. Wang, 64-75 

Zavala, J., Springel, V., & Boylan-Kolchin, M. 2010, MNRAS, 

405, 593 
Zhang, L., Weniger, C, Maccione, L., Redondo, J., & Sigl, G. 

2010, Journal of Cosmology and Astro-Particle Physics, 6, 27 



A. J. Cuesta et al. 



fc^ 

.s^* 



yiyi^cSirim 



Filament 1 // 









Fig. 3. — S/N all-sky maps from Fermi simulations for DM 7-rays in the energy range 100 MeV-10 GeV built from the Boxl60CR 
constrained simulation of the Local Universe. Results for DM decay (left) and annihilation (right) are shown for the bb channel model. 
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TABLE 2 

PHOTON COUNTS AND BACKGROUNDS IN Fermi SIMULATIONS. 



Object 


b b channel 
ann dec 


ann 


channel 
dec 


background 


Coma 1 deg 
Coma 2 deg 
Coma 5 deg 




5.297 
5.068 
3.245 


(44) 
(68) 
(85) 


5.297 

7.583 

10.078 


(44) 
(114) 
(303) 


1.237 
0.818 
0.445 


(7) 

(9) 

(11) 


2.109 
2.741 
3.245 


(13) 
(33) 
(85) 


25 
112 
601 


Virgo 1 deg 
Virgo 2 deg 
Virgo 5 deg 




5.041 
5.831 

4.185 


(36) 

(68) 

(102) 


5.646 
10.025 
15.588 


(43) 
(147) 
(488) 


0.000 
1.132 
0.666 


(0) 
(10) 
(15) 


1.118 
3.116 
5.068 


(5) 

(31) 

(126) 


15 

68 

492 


Perseus 1 deg 
Perseus 2 deg 
Perseus 5 deg 




0.777 
0.503 
0.271 


(10) 
(14) 


4.178 
4.042 
3.156 


(51) 

(88) 

(168) 


0.298 
0.152 
0.077 


(3) 
(3) 
(4) 


1.323 
1.041 
0.807 


(14) 
(21) 
(42) 


98 

386 

2665 


GAttractor 1 deg 
GAttractor 2 deg 
GAttractor 5 deg 


0.175 
0.130 
0.089 


(5) 

(8) 

(15) 


1.741 
1.654 
1.379 


(51) 
(103) 
(234) 


0.070 
0.033 
0.024 


(2) 
(2) 
(4) 


0.592 
0.438 
0.343 


(17) 
(27) 
(58) 


807 

3777 

28572 


Filament 1, d = 


= 65Mpc/h 


0.224 


(14) 


4.485 


(290) 


0.112 


(7) 


1.379 


(87) 


3891 


Filament2, d = 


: 40Mpc//i 


0.517 


(46) 


6.541 


(602) 


0.135 


(12) 


1.797 


(161) 


7869 


Filaments, d = 


= 65Mpc/h 


0.226 


(54) 


3.117 


(750) 


0.079 


(19) 


0.756 


(181) 


57127 


Filament4, d = 


= 55Mpc/fc 


0.380 


(60) 


5.486 


(881) 


0.120 


(19) 


1.338 


(212) 


24904 


Superclusterl, 


d = 45Mpc//i 


0.640 


(101) 


8.915 


(1445) 


0.177 


(28) 


2.343 


(372) 


24829 



Note. — The S/N ratio and number of photon counts (in brackets) in the 1 GeV— 10 GeV energy range 
for our different DM models. For cluster regions, three different radii are considered (1, 2, and 5 degrees). 
Filaments 1 to 4 represent elongated regions connected to these clusters which are potentially interesting 
due to their high S/N. Median distance of halos belonging to these filaments is indicated. Superclusterl is a 
collection of massive halos which accidentally lie along the line-of-sight. Background counts from the Galactic 
plus extragalactic diffuse in the same regions are also listed. Note that the annihilation to 6 6 case is shown 
for comparison purposes only. 

6. APPENDIX: ERRATUM TO THE PAPER 

The authors corrected a mistake in the computation of the area assigned to each pixel in the signal maps, which was 
missing a factor of 0.57rcos(6). While this error is not very important at mid-latitudes, the area of the pixels at the 
Galactic equator is 57% larger than previously estimated (hence lowering the signal), and is a factor of 2.5 smaller at 
the position of Virgo (6 ~ 75°), increasing the previously estimated signal arbitrarily as we approach to the Galactic 
poles. The signal-to-noise values of Table 1, as well as the maps in Figure 2 of the original paper are affected, and we 
report above the corrected version of both the Table and the Figure. Note that the dark matter density and density 
squared maps of Figure 1 of the original paper arc not affected by this error. The conclusions in the original paper 
remain valid and the values of the signal-to-noise for the high-latitude objects we focused in the paper (mainly Virgo 
and Coma) arc now even stronger. 

We have made the updated FITS files for the dark matter densit y and density squared maps available online at: 
http: //www, clues-project ■org/axticles/daj:kmattermaps.htinl[ 



AJC would like to thank the corresponding authors of http: //arxiv.org/abs/1110. 3305 (Gomez- Vargas et al., 
submitted to JCAP) for realizing this error when doing the analysis of their signal-to-noise maps. 



